基于多尺度方法的平动圆柱贮箱航天器刚——液耦合动力学研究

以充液航天器为工程背景,借助多尺度方法研究刚--液耦合动力学系统非线性动力学特性.利用多维模态方法,将描述横向外激励下圆柱贮箱中液体非线性晃动的自由边界问题转换为液体模态系数相互耦合的有限维非线性常微分方程组.推导液体晃动产生的作用于贮箱壁的晃动力和晃动力矩的解析表达式,进而建立航天器刚体部分平动和液体晃动耦合的非线性动力学方程组.应用多尺度方法对刚--液耦合系统的动力学特性进行解析分析,通过固有频率的特征方程求解耦合系统固有频率,推导外激励频率接近耦合系统第一阶固有频率时液体晃动稳态解的幅值频率响应方程.结合数值方法,研究了液体晃动稳态解的幅值频率响应曲线和激励--幅值响应曲线.结果表明, 随充液比变化,液体晃动稳态解的幅值频率响应曲线会发生软、硬弹簧特性转换现象和"跳跃"现象;幅值频率响应曲线的软、硬弹簧特性转换点受重力加速度和弹簧刚度系数影响;以上所得研究结果表明,考虑非线性效应时的刚--液耦合系统动力学特性与传统的线性系统模型所显示的动力学特性具有本质区别.本文的研究工作对进一步分析充液航天器刚--液耦合非线性动力学特性具有重要参考价值.      

接地式三要素型动力吸振器性能分析

利用固定点理论优化接地类型的动力吸振器得到的结果可能不是全局最优参数,在选择其他参数时主系统可以获得更小的振幅, 接地类型动力吸振器的优化问题值得进一步研究. 因此,以一种接地式三要素型动力吸振器为对象,通过研究系统参数变化对固 定点位置与主系统最大振幅的影响,得到了此吸振器的局部最优参数并分析了它的性能. 首先建立了此系统模型的运动微分方程, 得到了主系统振幅放大因子,发现系统存在3个与阻尼无关的固定点. 固定点中幅值较大点随系统参数变化的趋势可以代表最大振 幅随系统参数变化的趋势,因此利用盛金公式得到了固定点幅值的表达式. 为了更加精确,进一步使用数值算法得到了最大振幅与 系统参数的关系图,发现系统中存在局部最优参数. 通过对比接地式吸振器与接地三要素型吸振器的最大振幅随系统参数变化的趋 势,得到了接地式三要素型吸振器的局部最优参数,并发现当固有频率比小于局部最优频率比时,接地式三要素型吸振器模型主系 统的最大振幅要远小于接地式动力吸振器模型.      

主动射流控制水翼空化的数值模拟与分析

为了改善高速流动工况下水翼吸力面上流场的空化特性,提出了水翼表面主动射流对绕水翼周围流动加以控制的方法.基于密度分域滤波的FBDCM混合湍流模型联合Zwart-Gerber-Belamri空化模型,分析了来流空化数为0.83,来流攻角为8°,射流位置距水翼前缘为x=0.19c时,主动射流对于水翼吸力面上流动的空化特性和水动力特性影响.对回射流的强度进行了量化分析,以探究回射流与流场空化特性的关系.数值分析结果表明,在射流水翼吸力面上的时均空泡体积为原始水翼的1/15,使得流场内空化流动由云空化状态转变为较为稳定的片空化状态,显著地削弱了云空化的发展.此外,射流极大地改善了水翼的水动力性能,使得水翼的升阻比较原始水翼提高了22.9%,空泡的脱落频率减少了26.2%,空泡脱落所引起的振幅减小了9.1%.射流大幅降低了水翼吸力面上低压区面积,水翼吸力面上流体的逆向压力减小,回射流强度降低;同时,射流使水翼吸力面上的边界层减薄,增强了流动的抗逆压梯度能力,一定程度上阻挡了回射流向水翼前缘的流动,这也从机理上分析了主动射流抑制空化的原因.      

具有恒定冲击载荷的梯度泡沫金属材料设计

多胞材料可通过大变形大量地吸收冲击能量，引入密度梯度可进一步提高其耐撞性。梯度多胞材料的宏观力学响应对材料密度分布极为敏感，不同类型的细观构型的影响也极为不同。已有的研究工作主要局限在对给定的密度梯度分析其动态响应，较少对耐撞性设计方法进行研究。本文针对梯度闭孔泡沫金属材料，基于非线性塑性冲击波模型发展了耐撞性反向设计方法，以维持冲击物受载恒定为目标，运用级数法获得了简化模型和渐近解。利用变胞元尺寸法构建了连续梯度变化的三维Voronoi细观有限元模型，并利用ABAQUS/Explicit有限元软件对理论设计进行数值验证。结果表明，反向设计理论简化模型的渐近解对于梯度闭孔泡沫金属材料的耐撞性设计是有效的，所提出的耐撞性设计方法在控制冲击吸能过程和冲击物受载方面具有指导意义。     

基于双状态Χ~2检测和快速强跟踪AUKF的组合导航算法

针对GNSS/SINS/摄影测量组合导航中某个子系统发生故障时,整个导航系统易受到故障数据污染的问题,提出了一种基于快速强跟踪AUKF的双状态卡方(Χ~2)检测数据融合方法。首先,采用快速强跟踪AUKF算法进行滤波;然后,引入卡方检验通过检测UKF子滤波器输出的状态向量来定位故障参数;最后,采用强跟踪滤波准确跟踪状态矢量突变以增强系统鲁棒性,并根据自适应因子实时调整预测协方差阵以修正增益矩阵,使滤波结果不受异常信息的干扰。将提出的改进算法与常规算法分别应用于无人机着陆导航系统,结果显示:与传统UKF相比,提出的算法得到的位置误差减少了62.6%以上;与强跟踪UKF相比,导航误差也至少减小了32.6%。     

基于双轴转台的捷联惯性导航系统8位置系统级标定方法(英文)

惯性测量单元(IMU)是捷联惯性导航系统的核心设备,IMU的精度直接影响捷联惯性导航系统的精度。标定精度和速度是评价IMU标定方法的主要指标。传统标定方法是先进行温度补偿,再利用高精度三轴转台进行位置和速率实验,共需耗费约4天时间。虽然能够保证IMU的标定精度,但该方法复杂耗时。为改善复杂耗时的传统标定方法并保证标定精度,以便于批量生产,提出了基于双轴转台的8位置系统级标定方法。该方法基于27维状态向量的卡尔曼滤波器,以导航速度误差为观测量。该方法基于带温箱双轴转台设计,每个温度点需40 min进行标定。仿真和试验结果表明,该方法能够标定21个IMU误差参数,且该标定方法与传统标定方法相比标定时间更短,导航精度更高。     

一种考虑界面相物性影响的界面单元法

界面是由复杂的界面相简化而成的,界面破坏实际是界面相材料的破坏。数值计算为了方便,如经典模型和内聚力模型等,都把很薄的界面相作无厚度化处理。导致只能考虑界面的面力,而无法考虑界面相内的应力(平行于界面方向的应力)。使界面失效准则先天性地排除了界面相内部应力的影响,从界面相材料失效机理的角度来看这是不够严谨的。本文将界面相材料等效为一种弹性连续体,由界面本构关系推导得到了一种新的界面单元。该单元具有界面参数易确定、对界面相物性可以进行等效描述等优点。通过商用有限元软件ABAQUS和用户子程序UEL实现了数值分析,并与直接物理模型的数值模拟结果进行对比,证明了本方法的简便及准确性。通过对不同界面相厚度结构的进一步分析,探讨了本文方法的可行范围。     

基于Hopf-Cole变换的Burgers方程初边值问题的Sinc-Galerkin法

利用Sinc-Galerkin法数值求解Burgers方程的初边值问题。首先,用Hopf-Cole变换将二阶非线性的Burgers方程变换为二阶线性方程,同时把第一类边界条件变为第二类边界条件。时间上的导数采用θ加权格式离散,空间导数采用Sinc-Galerkin法离散,端点处分别引入权函数处理变换后的第二类边界条件。最后,通过数值算例验证了Sinc-Galerkin法的指数收敛性,与精确解相比,本文构造的数值格式精度高,能够有效捕捉激波等物理现象。     

非Hermitian正定线性方程组的外推的HSS迭代方法

为了高效地求解大型稀疏非Hermitian正定线性方程组，在白中治、Golub和Ng提出的Hermitian和反Hermitian分裂（HSS）迭代法的基础上，通过引入新的参数并结合迭代法的松弛技术，对HSS迭代方法进行加速，提出了一种新的外推的HSS迭代方法（EHSS），并研究了该方法的收敛性.数值例子表明：通过参数值的选择，新方法比HSS方法具有更快的收敛速度和更少的迭代次数，选择了合适的参数值后，可以提高HSS方法的收敛效率.     

三次Gauss和及二项指数和的混合均值

我们用三角和的性质研究一类三次Gauss和与两项指数和混合均值的计算问题,并给出一个精确的计算公式.